Characteristic Relations of Type-I Intermittency in the Presence of Noise 
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Near the point of tangent bifurcation, the scaling properties of the laminar length of type-I in- 
termittency are investigated in the presence of noise. Based on analytic and numerical studies, we 
show that the scaling relation of the laminar length is dramatically deformed from for e > 

to exp{i|e| 3/2 } for e < as e passes the bifurcation point (e = 0). The results explain why two 
coupled Rossler oscillators exhibit deformation of the scaling relation of the synchronous length in 
the nearly synchronous regime. 
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Intermittency is the occurrence of a randomly alter- 
nating signal between long regular (laminar) phases and 
relatively short irregular bursts Q. It is considered to 
be important as one of the routes to chaos in nonlinear 
dynamics. There have been extensive studies to manifest 
the route in terms of experiment as well as theory 
The scaling properties of the laminar length were studied 
for the first time by Pomeau and Manneville in the Lorenz 
model j^]. Based on the renormalization group equation 
[RGE] some other authors also investigated them || 
Recently, it was reported that the reinjection mechanism 
is another important factor that dictates the scaling re- 
lation of the laminar length |7]|| . 

Because noise is not avoidable in real environments, 
consideration of it is important to studying the realistic 
properties of nonlinear dynamical system |)|. The char- 
acteristics of nonlinear dynamical system in the presence 
of noise have been investigated by several authors based 
on the Fokker-Planck Equation [FPE] and RGE 

[^],[| , since the studies of Brownian motion initiated the 
stochastic modeling of the natural phenomena || . 

The system without noise converges to fixed points 
when bifurcation occurs but it does not exhibit infinite 
laminar phase under the succeeding random perturba- 
tion. So there is the possibility that the scaling proper- 
ties show quite different features from those of the con- 
ventional ones. In this respect, the recent investigation 
fiof comes into our notice which observed that the scal- 
ing properties of the laminar length are deformed in the 
nearly synchronous regime of two coupled Rossler oscil- 
lators |11|. 

In this Letter, we investigate type-I intermittency in 
the presence of noise before and after the tangent bifur- 
cation. Solving the Fokker-Planck equation [FPE] (|,||,[|, 
we derive the scaling relations of the laminar length in 
the closed channel (i.e., e < 0) and present the results by 
using numerical solution and simulation. Based on them, 
we also explain the scaling relations that grow exponen- 
tially in the nearly synchronous regime in two coupled 
Rossler oscillators and eyelet intermittency [ p)||l5| ]. 

The local Poincare map of type-I intermittency in the 



presence of nois£ is described as the following difference 
equation 



X n +1 



(1) 



where a is the positive arbitrary constant, e the channel 
width between diagonal and map, and D the dispersion 
of Gaussian noise £„. In the long laminar region, we 
can approximate the difference equation to the stochastic 
differential equation as follows Ji]: 



x = -V'(x) + V2D£(t), 



(2) 
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where dot and prime denote the differentiation with re- 
spect to t and x, respectively, £(i) is the Gaussian white 
noise such that (£(i')£(*)> = 5(t' - t) and (£(*)) = 
and V{x) is the potential given by V(x) = — i 
where c is the integration constant. The above equation 
can be considered as the equation of motion of the point 
particle under the potential V(x) and random perturba- 
tion £(t). The relation between return map and potential 
is given in Fig. 1. In this figure the stable and the un- 
stable fixed points correspond to the extremal points of 
the potential. 



Return 

Map i"i \ -y^— / 




V(x) 



FIG. 1. The relation between fixed and extremal points for 
e < 0. s and u axe the stable and unstable fixed points, xi 
and £1X6 the lower and upper bounds of laminar phase, 
respectively and x~ and x+ are the extremal points of the 
potential. 
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From the above stochastic differential equation, we 
obtain the backward FPE by following the well- 

established procedure M as follows: 



dG(x,t) = _ yl dG(x,t) 



dt 



dx 



D 



d 2 G{x,t) 
dx 2 



(3) 



where G{x,t) is the probability density of particle at 
{x,t}. We obtain a mean first-passage time [MFPT] 
equation after integrating the above FPE with respect 
to time as follows 



1 = -V'(x) 



dT 

dx 



D 



d 2 T 
dx 2 ' 



(4) 



where T{x) is the mean escaping time defined by T(x) = 
(t) = — J °° t dGi Q t ,t%> dt under the boundary conditions 
that G(x, 0) = 1 and lim^oo G{x, t) = 0. The MFPT 
function T(x) is the average transition time from the 
reinjection to the escaping point of the particle under 
the potential V[x) and random perturbation. 

The general solution of Eq. (3) can be derived as fol- 
lows: 



that -jjj 



T(x) 



I dx' exp 
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dx' / dx" 



exp 



D 



(V(x')-V(x"))\, (5) 



where c is the integration constant, xi is the lower 
bound of the laminar phase, and x is the destination 
point of the transition. We can easily verify that Eq. (5) 
is the general solution for the MFPT equation by insert- 
ing Eq. (5) into Eq. (4). 

If noise is small enough such that D <C 1, the first 
term in the above equation is suppressed by the factor of 
1/ D and the second term becomes dominant. The second 
term is not integrable analytically. Then we can expand 
the potential at the extremal point x± approximately (see 



v"0c±). 



Fig. 1) such that V(x) « V(x±) 
0((x-x ± f). 

In that case, the MFPT function T(x) can be approx- 
imated as follows: 



The above solution is consistent with the formal one 
which was derived in the previous investigation by 
the FPE and RGE analysis §-§,§ such that (I) ~ 
e" 1 / 2 f(<J 2 /e 3 / 2 ). We remark the fact that there has been 
no explicit derivation in analytic form like Eq. (7). The 
analytic solution is important in analyzing various inter- 
mittent phenomena quantitatively. In particular we are 
interested in the mysterious deformation of the type-I 
scaling near phase synchronization regime [jl0| and even- 
tually show that it enables the theoretical understanding 
of those phenomena. 

After taking the logarithm on Eq. (7), we obtain the 
equation such that InT ~ — | In |e| + |e| 3 / 2 . Our main 
interests are in the far region (|e| 3> 0) from the bifurca- 
tion point because the transition from the intermittency 
to stable orbit occurs here. In this region the scaling is 
dominated by the second term such that InT ~ |e| 3 / 2 
(note that this exponential saturation rapidly forms from 
|e| > 1.0 x 1CP 3 because the noise is small enough such 
0.5 x 10 +6 in our simulation (see Fig. 3 (b)) ). 
The reinjection probability P(xi n ) was reported to be 
another important factor which affects the scaling rela- 
tions of the laminar length [||. And it was generally 
considered to obtain average laminar length . But in 
this investigation we only consider the fixed reinjection 
probability P(xi n ) — 5(xi n — A) to study the intrinsic 
scaling property of the system (note in all of our simula- 
tion we set the reinjection point A = x{). 

Note that the approximation procedure used in Eq. (6) 
is not applicable for e > and a transient region because 
the equation is a solution for e < where the potential 
has extremal points (see Fig. 1). At the far outside from 
the bifurcation point (i.e., e ^> D > 0), we take the limit 
D — > in Eq. (4) and obtain the conventional scaling 
relation of type-I intermittency T ~ -1= . For inter- 
mediate range e ~ D, we present the numerical results 
of the scaling relation of the laminar length in Fig. 3 (c) 
and (d). 

In Eq. (7), if we take D -> the MFPT function T -> 
oo. So the particle trapped in the well does not escape 
from it if the random perturbation is turned off. 



T(x) 



1 

D 



exp 



{V{x+)-V(x-))}x 



£dx' jT ^"exp j^L [V"(x + )(x' -x + f 
-V"(x-)(x"-x.) 2 ]}. (6) 



The extremal points are given by x± = ±\J —e/a in Eq. 
(2). In the far outside of the laminar phase (i.e., at the 
limit x — > oo and xi -C x_), we can perform the integra- 
tion of the quadratic exponent || and then obtain the 
following approximated solution of the MFPT equation: 
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|T| = 




\3L>Va 



J3/2 



for e < 0. (7) 
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FIG. 2. The numerical solutions of the MFPT equation, 
(a) the first differentiation of the MFPT function and (b) 
the MFPT function when a = 0.1 and D = 2.0 x 10~ 6 . The 
dotted, dashed, and long dashed lines are for e = —1.0 x 10 -5 , 
e = -2.0 x 10~ 5 , and e = -3.0 x 10~ 5 , respectively. 

The scaling relation of the laminar length can be ver- 
ified in numerical solutions and simulations. In Fig. 2, 
the numerically solved the MFPT function T(x) shows 
the typical kink shape ]T^| , thus we can know that T(x) is 
the good physical quantity reflecting the transition char- 
acteristics of intermittency from the laminar phase to 
chaotic burst. In that case we can define a topological 
index of the transition such that Q = |T(oo) — T(— oo)| 
jl4j]. The negative signature in Fig. 1 stems from the 
backwarding property of FPE of Eq. (3) and the ab- 
solute value of the MFPT function T{x) is the laminar 
length |,|. 

In the following presentation of the numerical results, 
we let xi — —1.0 and x u = 1.0 and the average laminar 
length is (I) = \T(x u ) - T(x t )\ (as given in Fig. 2 (b), 
T(x) rapidly converges to the constants outside the cen- 
ter, so that we can say that (/) is a kind of topological 
index as defined above). To confirm the scaling relation 
of Eq. (7), we not only perform a direct simulation with 
Eq. (1) but also solve the second order MFPT equation 
(Eq. (4)) numerically. 
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FIG. 3. The laminar scaling before and after the tangent 
bifurcation, (a) and (b) are for e < 0. (c) and (d) are for e > 0. 
The circles are simulation data from Eq. (1) and the solid 
lines are solution data from Eq. (4). The dashed lines in (b) 
and (d) show the slope 3/2 and -1/2 saturation, respectively 
and the dot-dashed line in (b) is an analytic solution curve 
(Eq. (7)) (the maximum and dispersion of the Gaussian noise 
are \£\ = 2.0 x 10" 3 and D = 2.0 x 10~ 6 , respectively). 



In Fig. 3, the circle points and solid lines are sim- 
ulation data from Eq. (1) and solution data from Eq. 
(4), respectively. The dashed lines in Fig. 3 (b) and (d) 
are 3/2 scaling expected from Eq. (7) and —1/2 scaling 
for e > 0, respectively. When we simulate Eq. (1), the 
approximated Gaussian noise is used (see the caption of 
Fig. 3 for details). The solution data agree well with the 
simulation ones. We note here that there are some shifts 
from the solution line when the uniform noise is applied, 
but the scaling behaviors are invariant in both cases. 

Fig. 3 (a) shows the scaling relation for e < 0. The 
figure, ln(Z) as a function of |e| 3 ^ 2 shows a straight line 
approximately to confirm (/) ~ exp{-p|e| 3 / 2 }. The expo- 
nent 3/2 appears more clearly when we obtain ln(ln(Z)) 
as a function of ln(e) as given in Fig. 3 (b). Thus we 
can verify the slope is eventually saturated to 3/2. This 
is the very scaling behavior obtained in Eq. (7) analyt- 
ically. The analytic solution is also given in Fig. 3 (b) 
as the dot-dashed line that is the plotting of Eq. (7) . In 
figures though we present the simulation data within the 
limit of numerical calculation of the map (Eq. (1)) the 
data well follows the deformation of the scaling behavior. 

In Fig. 3 (c), the conventional scaling (I) ~ holds 
for relatively wide channel region (e ^> D > 0) so that 
the slope -1/2 saturation can be obtained in Fig. 3 (d). 
As the channel width e becomes close to zero, the straight 
line begins to bend (Fig. 3 (c)) and after e passes zero 
point, the straight line reappears (Fig. 3 (a)). 

We now apply this analysis to two coupled Rossler os- 
cillators which are a good laboratory for the studying of 
type-I intermittency with random perturbation [p~0| p~2|| - 
It is important to discuss the correspondence between 
two coupled Rossler oscillators and type-I intermittency 
in the presence of noise to explore the origin of nearly 
synchronous phenomena. The two coupled Rossler oscil- 
lators are given as follows ||ToHl^l : 

±i,2 = -^1,22/1,2 - Zi,2 + e(x 2 ,i - x lj2 ), 

2/1,2 = wi )2 a;i,2 + 0.15yi, 2 , (8) 

Zl,2 = 0.2 + ^1.2(^1,2 - 10.0), 

where u>i t 2 — 1.0 ± 0.015. The phase difference between 
the two oscillators can be rewritten as follows: 



-(9) = F(9,e) + G(4> 1 ,h), 



(9) 



where, 



F{9,e)=u 1 



A 2 



sm t 



G{4>\, 4>2) — 0.15(siri(/>i cos 0i — sin<rj 2 cos fa) 

+(— sm0i - — sm0 2 ). 
A\ A 2 



where 9 



Ax, 



+ yf 2 and 0i,2 



arctan(j/i i 2/xi i 2). In the above equations we neglect the 
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fast fluctuation term which depends on e, in F(0, e) as 
discussed in Ref . [jlO| . 

As already discussed in Ref. |l(J , the potential of this 
system V(9) = — J d8F(9,e) shows the saddle node bi- 
furcation at e — e t (= 0.0276). If we identify G(0i, <fe) as 
the random perturbation term of Eq. (1), we can argue 
the scaling relation of the length of synchronization in 
nearly synchronous regime of these two coupled oscilla- 
tors are effectively similar to that of type-I intermittency 
with noise (Eq. (2)) (note: the dispersion of G{4>\,4>2) 
hardly depends on e for e t < e < e c ). Thus we remark the 
scaling of the laminar length is separated into two regions 
in the center of e = e t . So it has (I) ~ exp{|e t — e| 3 ^ 2 } 
for e > e t and (I) — |e t — e|~ 1//2 for e < e t like the previ- 
ously presented results (note: based on numerical studies, 
Ref. [[To) proposed the scaling (I) ~ exp{— \e c — el 1 / 2 } for 
e > e t (e c = 0.0286)). 

As in the above analysis, the persistency of the in- 
termittency is caused by the random perturbation £(t) 
in closed channel region. Thus it can be argued that the 
nontrivial 2ir phase jumps observed in Ref. jl0| originates 
from the random perturbation term G(4>i, fa), neglected 
in the discussion of Ref. , rather than from the term 
^ + 3^- We perform the numerical simulation for the 
scaling of the laminar length in both regions (e < et and 
e > e t ), and the results are presented in the Fig. 4. As 
we expected, the figures well agree with the previous the- 
oretical analysis of the scaling relation. 
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FIG. 4. The laminar scaling in the two coupled Rossler 
oscillators for (a) e > e t and (b) e < e t . 

In conclusion, the new scaling relation of type-I in- 
termittency is presented in the region after the tangent 
bifurcation. We observe that the conventional scaling re- 
lation (I) -1= of type-I intermittency holds only in rela- 
tively wide channel (e^> D > 0) and it begins to deform 
as e approaches zero. The scaling relation is eventually 
saturated by (I) ~ exp{i|e| 3 / 2 } (Eq. (7) and Fig 3. (a) 



and (b)) after the tangent bifurcation (e < 0). Such dra- 
matic deformation of the scaling relation stems from the 
persistency of the intermittency with the random per- 
turbation even though the system is in a state of closed 
channel region (i.e., e < 0). This is why the laminar 
length in closed channel region grows faster than that of 
the positive one fll2||i~3j ]. From these results, we can also 
explain why two coupled Rossler oscillators exhibit defor- 
mation of the scaling relation of the synchronous length 
in the nearly synchronous regime [ fl0| |. 
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